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A Hamiltonian approach is presented to study the two dimensional motion of damped electric 
charges in time dependent electromagnetic fields. The classical and the corresponding quantum me- 
chanical problems are solved for particular cases using canonical transformations applied to Hamil- 
tonians for a particle with variable mass. The Green's function is constructed and, from it, the 
motion of a Gaussian wave packet is studied in detail. 



I. INTRODUCTION 

The motion of particles in vacuum and diverse media 
with dissipation has been studied in classical and quan- 
tum physics since long time. An important class of such 
problems are those of free electric charge carriers in a 
material under external time dependent electromagnetic 
fields. Of particular interest is the dissipation of energy 
through the interaction of charged carriers with the lat- 
tice ions (phonons) of the material, the carrier to carrier 
interaction through a Coulombian potential and, eventu- 
ally, through radiation. 

In classical systems damping is often described by in- 
cluding a velocity dependent drag term in Newton's sec- 
ond law. However, the inclusion of dissipation phenom- 
ena in quantum mechanics requires special care since its 
building blocks, time independent Hamiltonians, lead to 
energy conservation. This shortcoming is remedied in the 
heat bath approach [l|, 0] by coupling the single parti- 
cle Hamiltonian with an infinite degrees of freedom sys- 
tem, e.g., an infinite collection of harmonic oscillators, 
to which the energy of the single particle is transferred. 
Even though the energy is conserved given that the single 
particle, heat bath and coupling Hamiltonians are time- 
independent it is difficult to handle calculations with the 
many degrees of freedom of the heat bathQ. The dy- 
namics of an open quantum system Q is often formulated 
in terms of a master equation for the density matrix, 
that allows to work only with the single particle degrees 
of freedom by adding extra terms to the Von Neumann 
equation. Notwithstanding, the change over time of the 
open quantum system, in general, can not be presented 
in terms of a unitary time evolution [3]. 

In this work, we treat the problem of energy dissipa- 
tion by means of a single charged particle time depen- 
dent Hamiltonian [5|-|7[ . In contrast to time independent 
Hamiltonians, in the time dependent ones the energy is 
no longer a conserved quantity and therefore they allow 
for the possibility of energy loss. In particular, we study 
a Hamiltonian with a variable mass term under a time 
dependent electromagnetic field. 

In classical mechanics friction is usually analyzed in- 
troducing an opposing velocity-proportional force. The 
equation of motion of the particle can be usually built 



without difficulties from Newton's second law of motion. 
For a one dimensional particle with mass m subject to a 
potential U one has 



dU 
dx 



= 0, 



(1) 



where x is the position of the particle and r is the collision 
time. 

A deeper dynamical analysis is reached when the 
Hamiltonian formalism is applied. In the special case 
of one dimensional movement described by Eq. ([1]) the 
dynamics of a particle may be expressed by the Kanai- 
Caldirola Hamiltonian 0-0] 



H = ^-e~ t/T + U(x)e t/T . 
2m 



(2) 



This Hamiltonian even allows for analytical treatment 
in some simple quantum mechanical systems as a free 
particle]!] (U = 0) and the harmonic oscillator [9l-fl4j|. 

A great deal of effort has been focused on the modeling 
of dissipation phenomena for a ch arg ed particle through 
time dependent Hamiltonians [TBI. fl6|. However, obtain- 
ing a Hamiltonian for a dissipative charged particle un- 
der electric and magnetic fields is not as straightforward 
as for the Kanai-Caldirola Hamiltonian @. The naive 
assumption of a damping force proportional to the veloc- 
ity does not lead to a Hamiltonian formulation, i.e., the 
Newton's equations of motion 



m 

mx H x + qBy — qE x — 0, 

r 
m 

my H y — qBx — qE y = 0, 



(3) 
(4) 



of a particle in perpendicular electric and magnetic fields 
E x \ 4- E y j and Bk respectively can not be obtained from 
a Hamiltonian approach. For this reason, from here on 
we call this the naive model. 

Nevertheless, as we shall demonstrate below, it is possi- 
ble to model dissipation by introducing a time-dependent 
mass in the Hamiltonian for a charged particle 

H = ~ ttt (p - qA) 2 + q<f> + V. (5) 



2m (t) 
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The aim of this work is to study the dynamics of a 
damped charged particle in the presence of time depen- 
dent perpendicular electric and magnetic fields by means 
of a time dependent Hamiltonian. We obtain the general 
solutions for the equations of motion for the classical, as 
well as for the quantum problem, via the reduction of the 
Hamiltonian to zero by means of a series of linear canon- 
ical transformations in the classical case and correspond- 
ing unitary transformations in the quantum mechanical 
one. Here it is important to stress that, in general, in a 
large kind of dynamical systems the number of constants 
of motions is not enough to reduce the Hamiltonian to 
zero [13] ■ In this work, it is assumed that the Hamilto- 
nian is at most quadratic in the canonical coordinates, 
so that A is at most linear in the generalized positions, 
but the scalar potentials can be quadratic. 

The well known classical and quantum dynamics for a 
constant or a variable mass charged-particle in constant 
perpendicular electric and magnetic fields are recovered 
from our analysis. 

This paper is organized as follows. In Sec. |TT]we review 
the role of time dependent masses in the Hamiltonian of 
charged particles interacting with electromagnetic fields. 
In Sec. IIIII we address the solution of the classical Hamil- 
tonian via canonical transformations. The quantum me- 
chanical problem is introduced in Sec. IIVI Unitary trans- 
formations are applied to reduce the quantum mechan- 
ical Hamiltonian in Subsec. IIV Al With the resulting 
time evolution unitary operator, the Green's function is 
derived in Subsec. IIVBI As an example we study the 
dynamics of a Gaussian wave packet under the action of 
the Hamiltonian solved in this paper in Subsec. IIV CI 
We conclude in Sec. |V]with a summary of the results. 



and leads to a kinetic energy given in terms of the mo- 
menta as 



1 



T=-(p-qA,g- 1 (p-qA)) 



= o£(0 l )« (i* - a- 4 *) fa - • ( 9 ) 



In this approach it is assumed then, that the media acts 
on the particle by means of an alteration of the metric 
corresponding to replace the constant mass of the particle 
by a time dependent effective mass. Only the flat diagonal 
case gij = 8ijm(t), with a time dependent mass, shall be 
studied in here. However, more general metrics could be 
introduced in this manner, for example to include space 
inhomogencitics 18], but they shall not be considered in 
this work. 

Let us here start with the classical Hamiltonian for a 
charged particle 



H=^-(p-qAf 



(10) 



with a time dependent mass m. The equations of motion 
obtained from (1101) are 



(11) 



dH \ _Pi__ qAj 
dpi J m ttl 

dx t ~ m 4- dx % (Pj q 3 

j 



-<1 



d<f> dV 
dxi dxi ' 



(12) 



written as the Newton's second law they take the follow- 
ing form 



II. HAMILTONIAN WITH A VARIABLE MASS. 



To study the above physical problems a geometric 
setting is adopted. Let the kinetic energy T be given 
by a smoothly varying family of Riemannian metrics 
{r,gr} — Ylij 9ij( x ->t)%i%ii parametrized by time t on 
a compact 2n-dimensional manifold. The Lagrangian is 
then 



L = T -V - 



qA ■ r, 



(6) 



in terms of the vector potential A, the scalar potential <fi, 
and a confining potential V. The Hamiltonian is given by 
the Legendre transformation of the generalized velocities, 

H = ^2piXi-L = T + V + q<f>-qA-r, (7) 

i 

dL 



P.j 



(8) 



— (mr) = q(E + r x B) - VV. 

dt 



(13) 



with B = Vxi and E = —<p — dtA. It must be empha- 
sised that this equation is obtained from a Hamiltonian 
variational principle. A common approach is to set the 
time dependence of the mass through 



(14) 



m(t) = m e a[t \ 

where mo is a constant and a is a function of time [TBI. [l9|. 
From here on m shall label the time dependent mass 
and mo the constant one. Replacing this time dependent 
mass in the equation of motion (I13|) we obtain 

m e Q (r + df) = q(-Vc/>- d t A + r x V x A) - VV, (15) 
and, in terms of the constant mass m-o, 

m r = q(E a + r x B a ) - a(m r + qA a ) - W a . (16) 
Here, the new potentials 

A a =e- a A, (17) 

4>« = z~ a 4>, (is) 

V a = e~ a V (19) 
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have been introduced, and their magnetic and electric 
fields are identified as usual as 



B a =V x A a , 



dA a 
dt ' 



(20) 
(21) 



Although the expressions (|20|) and ([21]) for B a and E a 
are the standard ones, notice that the equation of motion 
(H"6f presents the extra term aqA a that could be inter- 
preted by considering that the damping is proportional 
to the generalized momenta rather than the velocities. If 
a > it might correspond to damping contributions. For 
a uniform magnetic field, similar terms have been found 
to be necessary to obtain a consistent quantum mechan- 
ical analog [HI, . 

A peculiar characteristic of the equations of motion 
(Unj) is that they are gauge dependent. This is natural 
since the (gauge independent) original physical law (TIB")) 
must be transformed using the minimal coupling that it 
is explicitly gauge dependent. Therefore, to start from a 
phenomenological model with constant mass, without the 
term aqA a in Eq. (|16[) . is not equivalent to the simplest 
physical model for damping. 

In order to illustrate how to model dissipation through 
a time dependent mass let us consider a charged particle 
in uniform perpendicular magnetic and electric fields 



B = Bk, 

E = EA - 



Eyl 



(22) 
(23) 



Separating the two components of Eq. (|13|) we obtain 
the following equations of motion for the particle 



mx + rax — muj c y — qE x = 0, 
my + my + muj c x — qE y = 0, 

where the cyclotron frequency 

qB 

u c = , 

m 



(24) 
(25) 



(26) 



is, in general, time dependent. In stationary state (x 
y = 0), the solution for these equations is 



q E x + (qB/m)E y 
m l + {q 2 B 2 /m 2 ) ' 
q E y - (qB/m)E x 
m 1 + (q 2 B 2 /m 2 ) 



(27) 
(28) 



In order to test the time dependent mass model equa- 
tions, specially the ones that describe the stationary 
state, let us try two different time dependent mass mod- 
els. First we consider a Kanai-Caldirola-like mass [l5l.[l6j 



m e 



t/7 



(29) 



where, for example, mo and r may be related to the 
effective mass and collision time in a semiconductor with 
mobility /i e = nq 2 r/mQ and charge carrier density n. 
Dislike the naive model, in this case, the time dependent 



mass model yields vanishing velocity components even in 
the presence of an electric field. Well known results, as 
the magneto conductivity tensor in semiconductors |20| . 
are contradicted by this calculation. 

As a second example let us consider the following con- 
venient choice of the mass' time dependence 



in = nu, [ - + k ) . 



(30) 



which we shall call the linear time dependent mass model 
(LTDMM, for "short"). Here A; is a dimensionless pos- 
itive parameter. For this particular case the stationary 
state solution for the velocity components is 



qr E x + LU TEy 
mo 1 + (jJqT 2 
qr E y - uj tE x 
too 1 + uj 2 t 2 



(31) 
(32) 



Thus, our LTDMM approach and the naive model given 
by Eqs. Q and ((4|) yield the same non vanishing sta- 
tionary state solution even-though their transient states 
might be different. 

To provide with a numerical example we have chosen 
a charged particle, e.g., an electron, in a GaAs sample 
with mobility [i e = nq 2 r/mo — 148m 2 /Vs, that yields a 
collision time r = 56ps. The effective mass and charge 
will be set to too = 0.067m e and q = — e, respectively, 
with e the electric charge of the electron. The magnetic 
and electric fields are B — 40mT and E = 100F/toj. 
The initial position and velocity of the particle are set 
to the origin and to f = x(0)i+ y(0)j = 3.7Km/ sj, 
respectively. 

Fig. [T] shows a comparison between the parametric 
plots of r (t) = x (t) l + y (t) j for the naive model (red) 
and the LTDMM (blue). We observe that even-though 
both models present different trajectories for the tran- 
sient state in t — > oo they have the same overall behavior. 

In Fig. [2] we can see a parametric plot of the veloc- 
ity vector r (t) = x(t)i + y (t)j, for the LTDMM (blue 
dots) given by (|3"D|) and the naive model (red solid line). 
Surprisingly both models plots are clearly over the same 
curve. Nevertheless we can not say that both examples 
behave exactly the same since the naive model reaches 
the terminal velocity faster than the LTDMM. This is 
shown in Figs. [3] and 2] where we can observe x and 
y plots for both models. We appreciate that the naive 
model saturates after t = 1 ns meanwhile the LTDMM 
saturates after t = 2.5 ns. 

Since the LTDMM yields similar results as the naive 
one, and both reach the same stationary state, we shall 
use it through out the rest of the work for the numerical 
examples. Notwithstanding, all the calculations in next 
sections do not rely on a specific mass model. 
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FIG. 1. (color online). Trajectory of a charged particle for 
the naive model (red solid line) given by Eqs.© and @ and 
the linear time dependent mass model (LTDMM) denned in 
Eq. (P0)> (blue solid line) and governed by Eqs.l[2"I)l and 




2 3 

x(t)(km/s) 

FIG. 2. (color online). Parametric plot of the velocity com- 
ponents for a charged particle for the naive model (red solid 
line) and the LTDMM for the classical case (blue dots) and 
the center of a quantum mechanical Gaussian wave packet 
(black crosses) as calculated in Subsec. IIV CI 



III. THE CLASSICAL PROBLEM: CANONICAL 
TRANSFORMATIONS 

A possible procedure to solve analytically the previ- 
ously depicted problem is to perform a set of canonical 
transformations [2l| • Equivalently, also the proposal of a 




FIG. 3. (color online). Velocity components, x (blue solid 
line) and y (red solid line), as functions of time for the naive 
model. 




FIG. 4. (color online). Velocity components, x and y, as 
functions of time for the LTDMM. 



function of the canonical coordinates at most quadratic 
in the momenta, has been successful in similar problems 
[22| . We chose this approach for the classical problem 
in order to establish a connection between the canonical 
and the unitary quantum transformation. 

The reduction of the Hamiltonian ^ is accomplished 
by applying canonical transformations of a certain sub- 
group of the afhne group, namely, translations, dilata- 
tions, shears, and rotations in phase space £ = {x,pY , 



(33) 
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with time dependent vector /x, and time dependent non- 
singular symplectic matrix M. 

The study of a charged particle's motion under homo- 
geneous electric and magnetic time dependent fields, is 
of utmost importance in the experimental and theoret- 
ical analysis of solid state devices. The building block 
of any theory explai ning the integer and fractional quan- 
tum Hall effects [23J, [24| , Shuvnikov-de Haas oscillations 
[25j . microwave induced resistance oscillations (26j, Hall 
induced resistance oscillations, amongst others, is the 2D 
electron in crossed electromagnetic fields. Therefore, we 
shall consider a 2D charge particle under perpendicular 
magnetic field 



B = Bk, 

with a vector potential given by 

B - B ~ 
A= -J yl+ 2**' 

The in-plane electric field is 



D 



-x — E„ 



with a scalar potential 



= -E x x - E y y. 



(34) 



(35) 



(36) 



(37) 



Here B, E x and E y are functions only of time. 

For the sake of simplicity and without any loss of gener- 
ality, we have considered the simplest gauge transforma- 
tion to write down the scalar and the vector potentials. 

The resulting quadratic time dependent Hamiltonian 
for the mentioned fields is 



" = ~ (pl+p 2 y) + I (rmv 2 c + «) {x 2 + y 2 ) 

--7r{ x Py-yPx)-qE x x-qE y yi (38) 



1 
2~^i 



for a confining potential V = k (t) (x 2 +y 2 )/8, and a time 
depending cyclotron frequency (|26[) . Since z is cyclic, the 
momentum associated with it has been dropped but not 
forgotten. Here all coefficients are given smooth functions 
of time. 

Our aim now is to reduce the Hamiltonian (1381) to zero 
using canonical transformations (see for example [27|) of 
a certain sub-group of the affine group. The procedure 
can be summarized as follows: 1) The third term in H , 
corresponding to the coupling z component of the angu- 
lar momentum, can be eliminated by a rotation leaving 
the first two terms invariant. The result is a Hamiltonian 
for two uncoupled one dimensional harmonic oscillators 
with variable masses and frequencies, as those considered 
in the literature; 2) a time dependent translation is per- 
formed to eliminate the linear contributions leading to 
an harmonic oscillator Hamiltonian with time dependent 
coefficients; 3) and, finally a dilatation and two shears 



are applied to reduce the Hamiltonian to zero; hence, 
the final generalized momenta and positions are simulta- 
neously constants of the motion and the original initial 
conditions, (j>o,qa), of our problem. 

For the first step in our program we require the gener- 
ating function of a rotation R for a finite angle 9(t) given 
by 

Fi = q'Rtpx 

= xp Xl cos 9 + yp Xl sin 9 — xp yi sin 9 + yp Vl cos 9, (39) 

with the column vectors q = (x,yY and p\ = {p Xl ,p yi ) i 
being pi — Rp and qi = Rq the rotated coordinates. By 
means of this generating function we obtain the following 
transformation rules 



Hi 



dFx 
dp Xl 
dF\ 

dPyi 



x cos 9 + y sins', 
—x sin 9 + y cos 6 



and 



Px = -q^- = p Xl cos 9 - p yi sin 9, 



dFx 
dy 



p Xl sin a + p Vl cos t 



(40) 
(41) 

(42) 
(43) 



Hence, the first transformed Hamiltonian is 



Of 



2m 

+ i {mul + k) (xl + yl)- qE* Xl ~ qE* yi . (44) 

Here H R is the original function, but now expressed in 
terms of the transformed coordinates and the rotated 
electric field 



E x = E x cos 9 + E y sin t 



E, 



-E x s'm9 + E y cost 



(45) 
(46) 



In order to reduce the angular momentum term in Eq. 
(gU), we set 



Y 



(47) 



and we obtain the direct sum of two one dimensional 
harmonic oscillator-like Hamiltonians 

- qE* Xl - qE? Vl . (48) 

All linear terms can now be reduced via space and mo- 
mentum translations with the generating function 

F 2 = (xi - X x ) (p X2 - n x ) 

+ (yi - Kj) (P V2 - %) - S, (49) 
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that yields the following transformation rules 

x 2 = q = x 1 -X x , (50) 

V2 = j. — =yi-X y , (51 

PX! = = Px 2 - 7T X , (52) 

Pvi = = Pv* ~ %' ( 53 ) 
where aca, 2/2, Pk 2 an d 

the new variables and S is 
the action. Here, X x and X y are time dependent parame- 
ters for the translation in coordinates, meanwhile ir x and 
TT y are the corresponding ones for the momentum space. 
After transforming via F2 the resulting Hamiltonian is 



H 2 



— ( 2 



+ g { muj l + K ) ( x 2 + vl) ~ + A z ) /' • ■ 



V m 



(mWj + /t) X x — qE x - ir x 



X2 



1 



1)2 



— (tt 2 + 7T 2 ) + - (mw 2 + k) (A 2 + A 2 ) - qE x X x - qE*X y + X x n x + X y n y - S. (54) 



— mui 

8 



In Hamiltonian if 2 the coefficients of X2, 2/2, Px2 and 
Pj,2 correspond to the Euler equations of the classical 
Lagrangian 

Ll = i ^ + ^ + I ^ + K) (A ' + ^ 

- qE*X x - qE^X y + X x ir x + X y ir y (55) 



for the translation parameters A^Ay, tt x and n y . In order 
for all the linear coefficients to vanish we require that this 
Lagrangian be the solution of the Euler equations for the 
translation parameters: 



d dLx 


dLx 


dt d~k x 


dir x 


d dLx 


dLx 


dt diiy 


dTTy 


d dL x 


dLx 


dtdX x 


^dK 


d 


dLx 


dtdX y 


dXy 



— + K = 0, 

to 



dLx _ fiiy_ 
dir„ V m 



mui 



rn 

2 _L 



tow 



IE* 



(56) 
(57) 
: 0, (58) 
0. (59) 



Additionally, to remove the Lagrangian part, Lx — S = 
must be fulfilled and consequently S can be associated 
with the time derivative of the corresponding action. 

The transformed Hamiltonian H2 is thus simplified 
into 



Ho 



1 

2m 



(*& + J&) + I ( muj c + K ) (£ + Vl) ■ ( 6 °) 



The harmonic oscillator coefficient tow 2 + k can be ex- 
pressed in terms of new parameters as 



TOW 2 + K = TO WQe 2 ^ a , 



(61) 



where e 2/3 = f 2 (t) + Koe a g(t)/mQU>Q contains the explicit 
time dependence of the given magnetic field B = Bof(t) 
and the confining potential k — Kog(t). 

In terms of this variables, the Hamiltonian is rewritten 
as [H] 



As a next step we consider a dilatation and two shears. 
The generating function for such a transformation is 



6 2 I 

cos 5 



(x 2 Px 



V2Py 3 



e 7 tan 5 
2A 

-^(*f+ % 2 ) (63) 



with time dependent functions 7, 5 and A. i 7 ^ produces 



7 



the following transformation rules Notice that even though the generating function F3 in Eq. 

j (|53"f has multiple divergences when 8 = (2n — l)n/2, its 

_ dF 3 _ ^ e 7 tanj corresponding canonical transformation rules (I64l) -(f67 l) 

op X3 cos 8 A have non. 

dF 3 e^ 7 e 7 tanJ 

y3 = ^r = c^ y2 -^^^' (65) 

<9F 3 ei 7 

-p^ 3 - A tan 8 x 2 , (66) 



8x2 cos <5 
<9F 3 _ e* 7 
9j/2 cos<5^ 



<9F 3 e^ 7 

Pj, 2 = — — = tPv3 ~ A tan 5 y%. (67) Under F3, the new transformed Hamiltonian is 



H, 



m 



cos 2 5 



mow5 e 
4A~ 



+ 



sin 5 ] — 5 + sin 8 cos <5- 



f 2 2 ^ 
7^ (P*a +P OT J 



,, , A >2 TOo^e 2 ' 3 2 \ • A 

' — sin 8 H 7 cos — — smo cos 8— 

m 4A / A 



Ae" 7 



(*3 + % 2 ) 



+ 



A moWpe 2 ' 3 



m 



4A 



• c r 7 • 2 t-A 

sm cos + — — sin 8 — 



(X3P X3 + VsPya) • (68) 



In order to obtain a null Hamiltonian we set the cocf- fication of the previous coupled equations 
ficients of (p 2 3 + p 2 3 ), {xl + yfj and {x 3 p X3 + y 3 Py 3 ) to 



zero. We, thus, obtain the following system of coupled 
differential equations for the transformation parameters 



i) , ' ( — cos 2 8 + 
m 



4A 



-8 + sin 8 cos 8—, 
A 



,2^2,9 



= e Q sm 2 5 H 7 cos 8 

\m 4 A 

—8 — sin (5 cos 8^r, 
A 

= e-N-— + ^|f!! , |sinJcos5 
m 4A 

sin 2 8- 



7 . o „A 



(69) 



(70) 



(71) 



8 = ^oe' 3 a cosh 77, 



7 = ^oe 



13 a sinh 77 (tan 8 — cot 5) 
sinh rj tan 5. 



(73) 

(74) 
(75) 



For practical purposes the solutions of these equations, 
in the most general case, can be obtained by numerical 
methods. Nevertheless, it is possible to extract informa- 
tion from (173")) - (1751) by grouping the last three equations 
in a single hyperbolic one 



5 1 

„2 



7 ~ $ - V 
b 2 



1, 



(76) 



here a = u>oe l3 ~ a /2 and b = woe^~ Q / sin2<5. If we use the 
r\ function as a parameter, we can rewrite the hyperbola 
with the parametric functions (1731) and 



The solutions to this differential equations cancel the 
whole Hamiltonian H 3 . In such a case x 3 , 2/3, p X3 and p y3 
are constant in time and, therefore, they are constants of 
the motion. To simplify the structure of the differential 
equations and their solutions we propose 



A = -mowoe^, 



(72) 



where r\ is a time depending function, yielding a simpli- 



j — $ ~ V = 



woe 



f) — a 



sin 28 



■ sinh 77. 



(77) 



For a given problem with no parabolic potential, n = 0, 
only one of the branches contains the physical solution. 
Each branch is associated with a given rotating direction 
of the charged particle. 

In particular, the vertices of the hyperbola correspond 
to the constant magnetic field case. If we set ourselves 
in one of the vertices 8 = u>^e^~ a /2 and by comparing 
with (|T3"|) we obtain that 77 = and, consequently, $ = 
and 7 = 0. This is indeed the case when the magnetic 



8 



field is a constant, i. e. /3 = and 7 = 0. In this manner 
we find that with an appropriate time dependent mass 
model and the initial condition 5(0) = we can integrate 
5 = woe~ a /2 and obtain 5, the only relevant parameter 
under the conditions described above. 

It is well-known that the time reversal symmetry is 
broken by a constant magnetic field, even though we 
have a frictionless problem, this symmetry breaking is 
the cause of the existence two vertices. More generally, 
for a problem where the magnetic field is a function of 
time, the solution is given by another region at the hy- 
perbola branch. In other words, the hyperbolic behavior 
of Eqs. (|73 ]) - ([75|) is a consequence of the magnetic field's 
time reversal asymmetry. 

The last transformation gives the solution to the ini- 



tial problem describing the motion of a charged particle 
under the influence of the potentials (|57|) and ([551 where 
the electric and magnetic fields are only time-dependent 
functions. Under the previous three transformations, we 
find that X3, 2/3, p X3 and p V3 are constants along the clas- 
sical orbit followed by the particle. In other words, X3, 
2/3, Px 3 and p V3 are the initial conditions and we shall re- 
name them as x 3 = x , y 3 = y , p X3 = p XQ andp M = p yo . 

It is also possible to figure out a single canonical trans- 
formation after adequately collecting all the above con- 
tributions into the following form 

£ = MCo + M, (78) 
where M is a symplectic matrix given by 



M 



a b 
c d 



e 2 cos 9 cos (5 

sin # cos (5 
-e~ 2 A cos 9 sin 6 
-e~i A sin sin <5 



— e 2 sin 9 cos 6 ^- cos 9 sin S — ^-sin#sin<5 

e~ 2 cos 9 cos S sin 9 sin S cos 9 sin 5 

e~ 2 A sin 9 sin <5 e 2 cos cos (5 — e"? sin#cos<5 

— e~ ' A cos 9 sin <5 e^sin#cos(5 e^cos#cos<5 



(79) 



and 





" cos 9 


— sin 9 





" 




X x 




M = 


sin 9 



cos 9 





cos 9 




— sin 9 




Xy 


■ (80) 




— 0— 


0— 


sin 9 


— cos 9 — 




- -*u - 


1 



As an example, we consider the simplest case when 
the magnetic field and the mass are constants, mean- 
while both the confining potential and the electric field 
are absent. In such a case S = — 9 = w /2t, A = tooWo/2 
and there is no dilatation, hence 7 = and 77 = 0. By 
using equation (|78l) all the position and the momentum 
variables can be expressed as a function of time and the 
initial conditions 



X 




y 




Px 




-Pv - 





+ coswoi) 
— \ sin wot 
— ■j mo wo sin w$t 
jm Wo (1 — cosw £) 



\ sin wo t 
\ (1 + cos w t) 
jinouio (cos wot — 1) 
— ^rnoWo sin w^t 



sin copt 

mom 
cosupt — 1 

\ (1 + COS Wot) 

— h sin wot 



1 — COS OJpt 

m wo 

sin u)gt 
mouio 

i sin Wot 
\ (1 + cosw £) 





" xo 




2/o 











(81) 



This last result is consistent with the solution obtained 
directly from the Hamilton equations of motion. The 
motion described in the previous equations is periodic, 
with period T = 2tt /w and w a is the Larmor frequency. 
The periodicity can be deduced from the behavior of the 
block matrices a and d in (|T9"|) since they become unit 
matrices for t = T, meanwhile b and c become zero. The 
charged particle is moving around a circular orbit in the 

plane xy with radium r — \jv\ Q + Py / m o UJ o- Physically, 
the trajectories of the particles are curved due to the 
Lorentz force, nevertheless, when the magnetic field B 
is small, the motion of the particles is almost linear (r 
grows). For larger values of Bo, the particle's motion is 
highly curved (r decreases). The last feature is given by 



the off-diagonal block matrices b and c. 

It is important to notice that in the Hamiltonian H 3 
we can set the two first coefficients to wo/2 instead of 
zero as in Eqs. (l69l) and (l70l) . meanwhile we keep the null 
equation (|TTj) . In this case we obtain a Kanai-Caldirola- 
like Hamiltonian, but the equations that must be satisfied 
in order to obtain a solution are much more complex. 



IV. THE QUANTUM PROBLEM: UNITARY 
TRANSFORMATIONS 

The classical calculations presented in the previous sec- 
tion allow to set a framework for a quantum mechanical 
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analog of (|38j) through the Schrodinger's equation 

H\il>{t))=pt\il>(t)), (82) 
where, the quantum mechanical Hamiltonian is given by 

(83) 



H 



1 

2m 



p-qA) 



hi 

8 



Here, p t is the energy operator, i.e., pt — > ifrdt and x, y, 
p x and p y are the space and momentum operators such 
that 



x\x,y) =x\x,y) , 
y\x,y) =y\x,y) , 

Px \Px,Py) =Px \Px,Py) , 
Py \Px,Py) =Px \Px,Py) , 



(84) 

(85) 
(86) 
(87) 



where \x,y) and \p x ,p y ) are the space and momentum 
eigenstates respectively. The space and momentum op- 
erators follow the usual commutation relations 



[xi,Pj] = ihSij, 
as well as the energy operator and time 
[Pt,t] = ih. 



(89) 



The physical electric and magnetic fields E and B, re- 
spectively, are obtained as usual from the scalar and vec- 
tor potentials cj> and A by the relations (|34|) and (j36j) as 
we discussed in Sec. |TTJ 

The integration of the quantum mechanical problem 
follows the same path as the classical problem. The re- 
duction of the Hamiltonian is now easily achieved by uni- 
tary transformations [291432T ] , each one associated to one 
of the three classical canonical transformations applied in 
Sec. Mil Each reduction step has the following structure 

UHU^U\ip(t)) 

= Up t U^U\yj (t)) H' W (i)) = p t W (t)) (90) 

with H' = UHW - U [p t , W] , and \ip' (t)) = U ' \ip (*)). 
The Floquet operator is thus given by 

H = H-Pt, (91) 

and the Schrodinger's equation takes the compact form 

H\ij(t))=(H-p t )\iP(t))=0. (92) 



Our aim now is to study the Hamiltonian in Eq. 
for the particular case analyzed in Sec. IHII of a magnetic 
and a perpendicular electric fields of Eqs. (|3"4"|) and (f3"r?)) . 
Such fields can be obtained from the potentials in (|35|) 
and (|37|) . The quantum mechanical potentials are thus 
given by 



A = — {-yi + xj) , 
4> = -E x (<) x - E y (t) y. 



(93) 
(94) 



In this gauge, the Floquet operator takes the following 
form 



n 



1 

2m 



Px 



qB 



Py 



qB 



q [E x (t) x + E y (t) y} + - (x 2 + f) - p t . (95) 



A. Evolution operator 

To obtain the evolution operator for the Hamiltonian 
in (|83| in the presence of the magnetic and electric fields 
given by Eqs. (|34|) and (|36|) . respectively, we proceed in a 
similar fashion to the classical case in Sec. Mil We apply 
a series of unitary transformations, each corresponding 
to a canonical transformation of the classical case. 

The first unitary transformation, a rotation around the 
z axis [29j], corresponds to the canonical transformation in 
Eq. (l39t and is given by 



Ui = exp 



.0L Z 



(96) 



where L z = xp y — yp x is the angular momentum along 
the z axis. It has the following effect on the position, 



momentum and energ 


y operators 




UixU\ 


= x cos 6 — ysin6>, 


(97) 


UiyUl 


= x sin + y cos 6, 


(98) 


u 1Px ul 


— p x cos 9 — Py sin 9, 


(99) 


UlPyU\ 


= p x sin 9 + p y cos 9, 


(100) 


UiptU\ 


= p t - 9L Z . 


(101) 



Note that U\ leaves invariant the quadratic forms x 2 + y 2 
and p 2 + Py, yielding the transformed Floquet operator 



UiHVt = ± (pi 



fy) 



1 



(E*x + EXy)-[8 



Pl . ( I D2 i 



where E x and Ey are the rotated components of the 
electric field given in Eqs. (|45l) and (14^1) . If the time 
dependent parameter 9 of the U\ transformation follows 
Eq. (|4"T1) it is possible to reduce the term proportional 
to the angular momentum L z . The Floquet operator is 
thus completely separated into the x and y parts. Now 
the Schrodinger equation takes the shape of two uncou- 
pled one dimensional harmonic oscillators. Here it is im- 
portant to set 9 (0) = as the initial condition for the 
parameter in order that U\ goes to unity as t — > 0. Wc 
will set this initial condition for all the transformations' 
parameters. Once the Floquet operator is separated we 
can proceed to reduce each part with the unitary trans- 
formations. We note that if m is time independent then 
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8 = —ujQt/2 with wo = qBo/mo, the cyclotron frequency. 
In this case, the charged particle motion is taken to a 
reference system that turns at half the cyclotron angular 
frequency. 

The next unitary transformation corresponds to dis- 
placements in space, momentum and energy and is asso- 
ciated with the canonical transformation (|49p . It is given 
by 



where 



Hit = exp 



U 2 = U 2t U 2x U 2y , 



-S(t) 



(103) 

(104) 
(105) 
(106) 



with time dependent transformation parameters S (t) , 
X x (i), (t), X y (t) and ir y (t). This unitary operator 
(|103p yields the following transformation rules 



U 2x = exp 


i 
n 


(t)x 


exp 


~zX x 
n 


(t)Px 


U 2y = exp 


i 


m 


exp 


"s 


(t)Py 



U 2 xU\ 


= x + X x (i) , 


(107) 


V 2 p x V\ 


= Px - K x (t) , 


(108) 


U 2 yUl 


= y + X v (£) , 


(109) 


U 2 p y u\ 


= Py ^x (t) , 


(110) 


v 2 p t v\ 


= Pt + S - X X TT X - XylTy 






+TT X X + X x Px + T^yV + XyPy. 


(111) 



Now we apply successively transformations Ui and U 2 to 
the Floquet operator (f9"Tj) obtaining 



U 2 UiUU\ul 



1 

2m 



(Px - n x ) 2 + (p y - ir y y 



+ i W + k) (x + X x ) 2 + - W + k) (y + X y ) 2 
o o 

- q [E% {x + X x ) + Ey (y + X y )] -p t -S (t) 

+ X x TTx + XylTy - TT X X - TT y y - X x px ~ Xypy. (112) 

In the previous transformed Floquet operator, as in the 
classical Hamiltonian, we identify the Lagrangian L\ of 
the transformation parameters and the corresponding 
Euler equations (|56"1) - (|59")l . Eq. (|112[) can be recast in 
the following form 



u,u 2 huM = ^(Pl + f y )+l W + K ) i* 2 + v 2 ) 



dtdXx 



dt dX„ 



dLi 

dX x 


x + 


' d dLi 

dt dn x 


dL{ 
dir x _ 


Px 


8Xy 


y + 


' d dL x 

dt dn y 


dLi 

dTTy_ 


Py 



Pt + L-S. (113) 



In order to reduce the linear terms and simplify the Flo- 
quet operator, we assume that the Euler equations (|56j) - 
([59)) are met for the parameters Aa;, X y , n x , ir y and S. 
The transformed Floquet operator (|112p is thus simpli- 
fied into 



Pi) 



k) 



Pi- 



rn) 



Corresponding to the F 3 canonical transformation in 
Eq. (|63p. the last unitary transformation can be split 
into the x and y parts as shown below 



U 3 = U 3x U 3 y. 



(115) 



The first unitary transformations in the right hand side 
is devoted to reducing the quadratic terms in the x part 
of the Floquet operator and correspondingly the second 
term reduce the y part of the Hamiltonian. 

The first transformation corresponds to a shear and is 
given by 



U 3x = exp 



4 H + ifi 



x exp 

and yields the following transformation rules 



+ 2 » e 2 

U 3x xU 3x = e 2 x cos 5 + ~^Px sin 6, 

U 3x p x ul x = e^px cos 6 — Ax sin 6, 
' A 



(116) 



(117) 



U 3xPtUl x = p t 



2A 



• sin 2 6 



(xpx + p x x) 




A 
2A 



sin 26 p\ 



Ae " 7 i A • or -2 
— [ S+ 2A Sm25 ) X ■ 



(118) 



In order to compute Eq. (|118[) it is necessary to ob- 
tain the time derivative of the unitary transformation. 
Although it can not be computed by direct derivation 
because the generator of (|116[) does not necessarily com- 
mute with its time derivative, it may be easily performed 
by separating the transformation (|116|) into a shear and 
a dilatation 



U: 



exp 



2h 

x exp 



(xp x + PxX) 
X / 

A Q x 2 



iS 
2H 



1 

A~ 



p: 



x exp 



m 

x exp 



(xpx+PxX) 



Ah 



(xp x + p x x) 



(119) 
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where A = Aoe 2/i and Ao is a constant. Here it is conve- 
nient to set the time dependence of the mass, magnetic 



field and confining potential by means of Eqs. (fT4"]) and 
([6~T]) . Applying this transformation to the Floquet oper- 
ator we readily obtain 



U 3x U 2 UiHUMul 



+ 



A 

to 



n c TO Wn e • 2 t- A • r c 

cos o H 7 sin o — o + — sin o cos o 

4A / A 



« ( — sin 2 6 
mo 



A TOw 2 e 2 ^ 



4A 



4A 

• a a A • 
sin o cos o — — sin 

A 



TO wge 2/3 2 \ • A 

■ cos o—o r sm o cos 



2 <5 



2A P * 
Ae -7 



i + + -^-pl + i (tow 2 + n)y 2 - p t . (120) 



r 



In order to vanish the terms proportional to xp x +p x x, p 2 
and x 2 , the differential equations for the 7, A and S pa- 
rameters between parenthesis should vanish. Notice that 
this equations are the same as Eqs. (|6"9")) - (fT0")) and conse- 
quently to Eqs. (|75| - ([75|) . Lastly the Floquet operator 
reduces to 

u 3x u 2 UiUuMul = -^fi 

+ ^(muj 2 c +K)f~p t . (121) 

The y part of the Hamiltonian can be eliminated by a 
transformation similar to Eq. (|116p given by 



U 3y = exp 



*2h 1 AlJ 



1 ~2 

A p y 



x exp 



-i-jftivPv +p y y) 



(122) 



In this transformation the parameters 7, A and S are the 
same as those from Eq. (|1 16[) since the x and y parts 
of the Hamiltonian are symmetrical. By applying this 
transformation we finally obtain 



it = U 3 U 2 U 1 HUlu^ul 



-Pt- 



(123) 



In Eq. (|123[) . the Floquet operator was reduced to the 
energy operator p t implying that any ket applied to the 
right of H as 

UU 3 U 2 XJ X |V (t)) = -ptU 3 V2Ui |V (t)) = 0, (124) 

should be a constant one, according with Schrddinger's 
equation ([92")) , e.g., 



/7 3 C/ 2 J7i |V (*)) = IV (0)) • 



(125) 



As a consequence, the state of the system at any time 
|V (t)) is connected to the state in t = 0, |V (0)), by 



|V(t)) = c/ 1 t c/Ml^(o))- 



The time evolution operator is thus easily obtained as 

U(t) = uMul (127) 

and the state of the system at any given time t evolves 
from t = according with 



|V(t))=w(*)|V(0)) 



(128) 



(126) 



Notice that time enters the evolution operator U through 
the parameters 9, X Xl n Xl X y , n y , 7, 5 and A in each of 
the unitary transformations. 

We can calculate the position and momentum opera- 
tors in the Heisenberg representation by performing the 
three transformations on Schrodinger representation of 
the space an momentum operators 

x H (t) = U^xU = U 3 U 2 UixU\ulul, (129) 

y H (i) = UHU = U^UnjUlulul (130) 

PxH (i) = U^p x U = U 3 U 2 U x p x Ululul, (131) 

p yH (i) = tfp v U = U 3 U 2 U!p y UlU^Ul (132) 

By working out the explicit form of the previous transfor- 
mations using the transformation rules (ji?T)) - (|101|) . (|107[) - 
(fTTT]) and (TIT7| - (fTT8)) we obtain 

xh (t) — e _2 cos 9 cos 5x — e - ' sinflcosoy 

e 2 _ e 2 

+ — cos 9 sin Sp x — — sin^sinopj,, (133) 

y H (t) — e~ i sin 9 cos Sx + e~ % cos cos Sy 

x 2. 
e 2 e 2 

+—sin9sva.8p x + — cos6sm.5p y , (134) 

PxH (t) = — e~ 2 Acos0sin&E + e~% A sin 9 sin 5y 

+e^ cos 9 cos Sp x — e 2 sin 8 cos Sp v , (135) 

PyH (t) = — e~ 2 A sin 9 sin Sx — e~ % A cos 9 sin Sy 

+ei sin 9 cos Sp x + e 2 cos9 cos Sp y . (136) 

Here, it is worthwhile noticing that, as in the classical 
case, the previous equations can be also expressed in the 
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symplectic form as 







X 


fjH 


= M 


V 


PxH 




Px 


-PyH - 




-Py - 



(137) 



where M and fj, are given by Eqs. (|79|) and (|80|) . respec- 
tively. 



The shear matrix element is in fact the propagator for an 
harmonic oscillator [33J 



27rfi,sin<5 



,- 2B i iji I 23131 ' - ( x2 + X ' 2 ) COS S ) 



(146) 



B. Green's function 

The Green's function is calculated as usual in terms of 
the evolution operator as 



G(x,y,t\x',y',0) = (x,y\U (t)\x',y'}. 



(138) 



To obtain the explicit form of G we first calculate the 
matrix elements of each of the unitary transformations 
U\, Ui and U3 in order to join them by the integral 



G (x, y, 1 1 x', y', 0) = / / dx\dx\ (x, y U{ xi,y^ 



x (xi,yi 



Ui 



X2,V2 {X 2 ,V2 



x',y' 



(139) 



For Ui and U2 it is convenient to explore their effect 
on a space eigenstate. The rotation U\ has the expected 
effect on any space eigenket 

U{ \x,y) — |xcos6* — ysiri9,y cos9 + x sin 9) , (140) 
and its matrix element is hence given by 

x, y\ — 5 (x' — x cos 9 + y sin 9) 

x 5 {y' — ycos9 — xsin9) . (141) 



x',y' 



U 



The transformation U2 is a translation in space and mo- 
mentum, therefore its effect on a space eigenstate is 

U\ \x, y) = e-^e" 1 ^ \x + X x , y + X y ) , (142) 
and its matrix element is thus given by 



x',y' 



U. 



x,y 



» e 



x6{x'-x-X x )6(y'-y-X y ). (143) 



The transformation U3 is the product of a dilatation 
and a shear. The dilatation has the following effect on a 
space eigenstate 



e 4 



xe 2 



(144) 



where the coefficient e - 1 is due to the rescaling of space 
and the consequent renormalization of the space eigenket. 
The matrix element of the dilatation is then given by 



-i^(xp x +p x x) 



= e 4 S 



xe 2 



(145) 



After reducing all the integrals in (| 1391) the explicit 
form for the Green's function is obtained as 



G(x,y,t\x',y',0) 



S(t) 



27r?isin<5 

x e" 2 ^ co ^ cos s ^ n -^^ ) 2 +(y cos x s ^ n -^y )^ 
x e -|(' r =«+ijf7 a: ' e_7/2 )( a:co s 9 +»siii9-A a: ) 

x e --i;( 7r !/ + ^#7l'' c_ "' /2 )(l' cose - a:sin6 '- A !/) 

xe Tcoti(-' 2 +v' 2 ). 



(147) 



This Green's function has indeed the correct shape pre- 
dicted by Schwinger and others pjil - l^ : it should be com- 
posed only of linear and quadratic terms of the space and 
momentum operators. 



C. Gaussian Wave Packet 

As an example, we now wish to study the evolution of a 
charged particle Gaussian wave packet under the action 
of constant and uniform crossed electric and magnetic 
fields. For the mass we select the LTDMM from Eq. 
(|30p . and we set the same parameters from Sec. [TT1 in 
order to prove Ehrenfest theorem. 

We start (in t = 0) with a Gaussian wave packet of the 
form 



if>(x,y,0) 



I f x 2 + y 2 



exp 



2a 2 



, .Pxox+Pyoy . n .„, 
x exp J 1 — - — ) , (148) 



where p x o, p y o and a are the initial momentum and wave 
packet width values. Note that initially the wave packet's 
center is located at the origin, and since the constant 
magnetic field is directed along the z axis, the vector 
potential is given by (|3"5)) therefore its average vanishes. 
In this manner, the initial momentum and velocity are 
related by p x0 i + p y o3 = mf (0). 

The wave function at any time is explicitly calculated 
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in terms of the transformation parameters as 



i> {x, y t t) = J dx'dy'G (z, y, t \ x', y', 0) i/> (a/, y', 0) 
e _ "s s ( t ) / fte2 sin (5 



aA 



+ iae 2 cos <5 



x exp 



1 

2^2 



x exp 



x exp < i 



2na 2 
. A cot 8 



(x-c?) +{y-0 



21i 



(x-X R Y + (y-Xl 



x exp -1 (x - A*) + < („ - A«)] , (149) 



where the standard deviation of the wave packet is given 
by 



a(t) = \ /a 2 e-T cos 2 S + 



h 2 e~> sin 2 5 
a 2 A 2 ' 



(150) 



and correctly complies with a (0) = a. The rotated A, 7r 
and £ parameters are given by 





= X x cos t 


) — Xy sin0, 


(151) 


K 


= Xy COS t 


' + X x sin#, 


(152) 


n x 


= n x cos 6 


' — Tt y sin#, 


(153) 


n y 


= TTy cos £ 


' + TT y sin#, 


(154) 


S.r 


= (x COS 6 


- Cy sin6>, 


(155) 




= C;/ cos 6 




(156) 



The C parameters are composed of two parts 

Cx = X X + Aj;0, 

Cy = A y + XyQ , 

where 



A x o — 



2 



AyO - 



A esc 6 
_ 

ea 

A esc 5 



PxO, 
PyO- 



(157) 
(158) 



(159) 
(160) 



The probability density can easily be worked out from 
Eq. (I149P giving 



\^(x,y,t)\ 2 = — ^exp 



' a 2 \ X Cr 



x exp 



(161) 



From the previous expression it is clear that the Gaus- 
sian wave packet follows the trajectory given by the 



vector C — ( R i + Moreover, using the differen- 

tial Eqs.flU}, t[72 ]) -([75 |) and (HHJ) - (ESJ) it is easily demon- 
strated that A^ and A^ fulfill the same equations of mo- 
tion as the classical particle 



iX 



'' mu; c X R 



mX R - qE x = 0, 



mXy + muj c X R + mAj - qE y = 0, 
and A^ and Xy fulfill the homogeneous equations 



(162) 
(163) 



mX 
mX 



R 



x0 
R 

yo 



mu c Xy - 
- mui c X R Q 



rhX R = 0, 
- mX R = 0. 



(164) 
(165) 



A y +A y0 



We can thus infer that ( R = X R + X R and C, R 
are the complete solutions for the classical equations of 
motion where X R and X R are the particular solutions of 
the inhomogeneous equations and X R and A^ are the 
homogeneous solutions baring the initial conditions. 

In this manner, the center of the wave packet follows 
the same trajectory as the classical particle. This is a 
proof of Ehrenfest theorem. The trajectories obtained for 
£ R indeed are the same as the classical ones as was proved 
by direct numerical calculations of the wave packet center 
motion shown in Fig. ((2|) with black crosses. 



V. CONCLUSIONS 

To summarize, we have studied the classical and quan- 
tum dissipation of a charged particle in variable mag- 
netic and electric fields through a time dependent mass 
Hamiltonian. To integrate the classical Hamiltonian, a 
series of three canonical transformations are explicitly 
constructed and applied in order to reduce it to zero. 
The final transformed variables are at the same time con- 
stants of the motion and initial conditions for the gener- 
alized momenta and positions. The final solution to the 
equations of motion is rendered in its symplectic form. 
Correspondingly, the quantum Hamiltonian is reduced 
to zero by three unitary transformations. This proce- 
dure allows for the calculation of the evolution operator 
in rather general conditions, i.e. time dependent mass, 
variable electric and magnetic fields. The generalized mo- 
mentum and space variables in the Heisenberg picture 
are expressed in terms of a symplectic linear combina- 
tion of their Schrodinger picture versions. In times, the 
Green's function is constructed from the evolution oper- 
ator and the calculated expression is consistent with the 
structure obtained by Schwinger and others [H|-[36j]. As 
an example, the dynamics of a Gaussian wave packet un- 
der damping and constant crossed electric and magnetic 
fields is studied. Its motion is proved to follow the same 
trajectory as the classical particle under the exact same 
conditions. The results presented in this paper might be 
useful in solid state calculations where dissipation plays 
an important role. 
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